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CHAPTER  I 
INTRODUCTION 

| 

Exploding  wires  with  high  current  pulses  have  been  studied  for 

about  two  decades J In  recent  years  pulsers  capable  of  delivering  on 

the  order  of  terrawatts  of  power  for  tens  of  nanoseconds  have  become 
3 

available.  These  high  power  pulsers  have  been  used  to  produce  electron 

4 5 

beams  and  to  explode  wires. 

The  Gamble  II  generator,  for  example,  consists  of  a Marx  generator 
which  stores  228  kJ,  an  intermediate  storage  capacitor,  a pulse  forming 

t 

line  and  a coaxial  impedance  transforming  line  with  an  output  impedance 
of  1.5  ohms.  It  has  been  used  to  deliver  on  the  order  of  40  kJ  to  an 
exploding  tungsten  wire.5’6 

Although  such  a device  produces  tens  of  Joules  of  radiation  in 
the  soft  X-ray  spectrum,  it  is  of  interest  to  examine  the  possibility 
of  scaling  the  technique  to  higher  power  levels.  The  usefulness  of  a 
scaled-up  power  supply  to  increase  the  soft  X-ray  radiation  depends 
upon  how  well  the  resulting  plasma  load  matches  the  high  power  pulser. 


Exploding  Mire  Plasma  Load  Characteristics 

Mosher,  et  al,  have  reported  that  "appropriate"  plasma  loads 
tend  to  match  to  the  generator  driving  them.5  This  was  found  to  be 
true  for  tungsten  wires  on  the  order  of  25y  diameter,  3 to  4 cm  in 
length,  open  circuit  power  supply  voltages,  vQ,  in  the  range  2-5  MV. 
It  remains  to  be  seen  whether  or  not  the  matching  will  extend  to 
appreciably  higher  power  levels. 
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The  power  supply,  or  pulser,  is  adequately  represented  by  a voltage 
source  and  a resistance  for  the  time  period  of  interest.  The  exploding 
wire  and  its  mounting  jig  can  be  represented,  insofar  as  its  effect 
on  the  circuit  is  concerned,  by  a time  varying  series  R-L  impedance.  The 
circuit  is  shown  in  Fig.  1. 


Fig.  1 Equivalent  circuit  for  exploding  wire  pulser 

The  simplicity  of  the  circuit  in  Fig.  1 is  highly  misleading, 
since  the  innocent  looking  terms  L(t)  and  R(t)  are  actually  extremely 
complicated  functions  of  the  exploding  wire  plasma  parameters.  The 
plasma  parameters,  in  turn,  are  sensitive  to  the  voltages  and  currents 
in  Fig.  1.  Because  of  the  complex  non-linearities  involved,  any 
theoretical  treatment  of  the  problem  must  involve  computer  aided  numer- 
ical solutions  for  specific  cases.  To  date,  no  computer  program  is 
available  which  describes  how  the  exploding  wire  plasma  reacts. 

5 

A one-dimensional  code  has  been  reported  which  includes  the 


radiation  and  electron  transport  equations.  Maxwells  equations  and  the 
external  circuit.  This  code  is  a significant  piece  of  work  which 


utilizes  the  Corona  equilibrium  model7  to  describe  the  local  high-temp- 
erature properties  of  the  plasma.  Although  the  results  obtained  with 
this  code  agree,  in  a broad  brush  sense,  with  the  experiments  performed 
to  date  with  power  levels  on  the  order  of  a terrawatt,  the  authors 
point  out  that  the  effect  of  run-away  electrons  has  not  been  included. 
And,  of  course,  since  it  is  one-dimensional,  it  does  not  include  the 
effects  of  the  observed  m = 0 instabilities. 

Apparently,  the  run-away  electrons  are  created  in  the  low  density 
regions  at  the  surface  of  the  plasma  cylinder.  Chapter  II  of  this  report 
considers  the  radiation  from  the  time  varying,  non-Maxwellian  velocity 
distribution  which  developes  in  a plasma  in  the  presence  of  a strong 
electric  field.  The  remainder  of  Chapter  I discusses  in  a semi  quanti- 
tative manner  the  load  R(t)  and  L(t)  characteristics,  the  m = 0 insta- 
bility and  the  scaling  up  problem. 

Kirchhoff's  voltage  equation  for  Fig.  1 may  be  written 

e(t)  = [Rg  + R ( t ) ] i(t)  + ^ [L(t)i (t ) ] (1) 

The  inductance  is  the  ratio  of  the  flux  within  the  plasma  plus  the 
flux  within  the  jig  but  outside  the  plasma  to  the  current.  This 
inductance  varies  with  time  for  two  reasons.  First,  the  flux  per  unit 
current  within  the  plasma  varies  because  the  current  density  distri- 
bution and  the  radius  of  the  plasma  both  vary.  Second,  the  flux  out- 
side the  plasma  per  unit  current  varies  because  the  plasma  radius 
varies.  (It  is  not  clear  from  the  reference  cited  above  that  this 
second  L variation  was  included  in  the  one-dimensional  code,  although 
it  may  have  been.)  It  will  be  seen  that  the  variation  due  to  the  second 
reason  are  larger  than  those  due  to  the  first. 
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Consider  the  inductance  due  to  the  flux  inside  a cylindrical  plasma 
with  a 2 mm  radius,  rQ,  and  4 cm  long,  with  the  current  uniformly 
distributed, 

L = (1/2  x 1 0-7 ) ( 0 . 04 ) = 2 x 10'9  hy  (2) 

= 2 nhy. 

Now  if  the  plasma  is  compressed  to  0.05  mm  radius,  rp,  the  inductance 
becomes 

p 

L = 2 x 10‘9  + 2 x 10"7  l In  ~ (3) 

P 

= 2 + 2(ln  40)  = (2+  7.38)  = 9.38  nhy 

If  the  current  density  within  the  plasma  is  higher  near  its  outer  radius 

and  lower  near  its  center,  the  inductance  in  eq.  (1)  and  the  first  term 
in  eq.  (2)  are  even  smaller  than  2 nhys. 

If  the  diameter  of  the  wire  were  25u,  the  inductance  due  to  the 
flux  inside  a 4 cm  cylinder  out  to  a radius  of  2 mm  is  12.15  nhy.  If 
the  wire  explodes  and  the  plasma  expands  to  a 2 mm  radius,  the  inductance 
drops  from  12.15  nhy  to  less  than  2 nhys  if  the  current  density  is 

larger  near  the  plasma  surface.  In  the  process  of  this  expansion,  some 

plasma  energy  is  removed  and  fed  into  the  magnetic  field.  If  the  plasma 
is  subsequently  compressed  to  a radius  of  0.05  mm,  the  inductance  goes 
back  up  to  a value  between  7.38  and  9.38  nhy,  depending  upon  the  current 
distribution.  During  the  compression,  some  magnetic  field  energy  is 
converted  into  plasma  energy. 

Of  course,  there  is  an  additional  inductance  due  to  the  flux  within 
the  jig,  but  outside  the  2 mm  radius,  ('50  nhy)  which  must  be  added  to 


I 
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the  above  inductances  to  obtain  the  L(t)  for  eq.  (1).  Numerical 
solutions  of  eq.  (1)  with  typical  R(t)  and  L(t)  values  and  with  e(t) 
a step  function  of  2 MV,  and  Rg  = In,  gives  a current  which  rises  to 
the  order  of  1 MA  in  a few  tens  of  nanoseconds.  In  some  cases  the 
current  does  not  monotonical ly  increase,  but  has  a slight  dip  in  the 
region  where  dL/dt  is  positive.  However,  since  the  current  is  higher 
during  the  compression  stage  than  during  the  expansion  stage,  a net 
energy  of  the  order  of  a kilojoule  can  be  fed  into  the  plasma. 

The  effective  resistance  of  the  plasma  is  very  important,  since 

2 

most  of  the  heating  is  due  to  the  I R loss  within  the  plasma.  Exami- 
nation of  the  approximate  resistance  of  a 4 cm  long,  25p  diameter 
tungsten  wire  as  a function  of  temperature  is  enlightening.  The  resistiv- 

O 

ity  of  tungsten  at  300°K  is  5.65  u ohm  cm  and  rises  to  117.1  at  3655°K. 
The  wire  resistance  goes  from  4.6  up  to  95  ohms  as  the  temperature 
increases  from  room  to  melting  temperature.  The  effective  resistance 
of  the  wire  load  may  not  reach  a value  as  high  as  95  ohms  since  some 
thermionic  electrons  will  be  emitted.  Also,  no  data  is  available  for 
the  resistance  of  the  tungsten  during  the  time  that  it  is  being 
vaporized  and  ionized.  Fortunately,  the  value  of  resistance  during 
this  time  frame  does  not  appear  to  be  too  critical  insofar  as  the 
high  temperature  matching  problem  is  concerned.  For  example,  if  the 
plasma  resistance,  R(t)  in  Fig.  1 were  10  ohms,  then  a large  fraction 
of  the  generator  voltage  would  be  dropped  across  R(t).  A voltage  of 
say  500  kV,  would  produce  a power  of  25  Gw,  which  would  represent  an 
energy  of  250  joules  in  10  ns.  Since  only  a few  tens  of  joules  are 
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needed  to  vaporize  and  ionize  the  tungsten  once  or  twice,  a resis- 
tance between  10  and  100  ohms  would  only  last  for  10  to  20  ns  at  most. 

5 

The  resistivity  of  the  ionized  plasma  above  about  10  °K 

g 

can  be  calculated  from  equations  from  Spitzer's  book.  In  the  absence 
of  a magnetic  field  the  plasma  resistiveity  is  given  by 


n 


3.8  x 103  Li&r  ohm 


J/2 


cm. 


(4) 


where  Z is  the  ionic  charge,  y£  is  a factor  given  in  Table  1,  T is  the 
temperature  in  °K,  and  ZnA  is  given  in  Table  2. 


Table  1 Factor  for  equation  (4) 
Ionic  Charge  1 2 4 16  °° 


ye 

0.582 

0.683  0.785 

0.923  1.000 

Table 

2 Values  of  ZnA 

for  equation  (4) 

Ion  density  (cm” 

3) 

Temp. 

<z> 

1019 

1020 

105 

6 

5.79 

4.64 

3.49 

!06 

11 

7.79 

6.64 

5.49 

107 

47 

9.37 

8.22 

7.06 

108 

63 

11.52 

10.37 

9.22 

The  average  values  of  Z in  Table  2 are  obtained  from  Mosher.'7 

Tn  a ctrnnn  mannptir  fiplri  trantvprcp  tn  tho  plprtrir  f-iolrl 
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the  resistivity  becomes 


nL  = 


1.29  x 10 


4 Z £n  A 

~^72~ 


ohm  - cm 


(5) 


The  calculations  below  indicate  that  a strong  transverse  magnetic 

field  exists  near  the  plasma  surface  at  high  temperatures  and  currents, 

but  not  in  the  central  regions  of  the  plasma. 

The  magnetic  induction  at  the  surface  of  a plasma  cylinder  of 

radius  r is  pQI/27rr,  which  yields  10  tesla  with  0.1  MA  at  r = 2 mm. 

Since  the  observed  radii  are  in  the  range  of  0.1  to  1 mm  for  currents 

on  the  order  of  1 MA,  there  is  a strong  magnetic  field  near  the  surface 

of  the  plasma.  The  resistivities  for  a tungsten  plasma  with  ion 
*18  20  -3 

densities  of  10  and  10  cm  with  a strong  transverse  magnetic  field 
are  given  in  Table  3. 


TABLE  3 Resistivity  of  Tungsten  Plasma  in  a 
Strong  magnetic  field,  (ohm  - cm) 


T 

00 

o 

II 

c 

n.  - ,020 

°K 

-3 

cm 

cm  3 

105 

1.42  x 10"2 

8.54  x 10'3 

106 

1.11  x 10'3 

7.79  x 10'4 

io7 

1.80  x 10~4 

1.35  x 10'4 

108 

9.36  x 10“6 

7.49  x 10'6 

The  resistance  of  the  plasma  would  be  given  by 


RD-G  = 


n 


l 

A 


(6) 


if  the  current  were  uniformly  distributed  over  the  plasma  cross  section. 
However,  the  resistivity  is  sufficiently  low  and  the  time  period  of  the 
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current  pulse  is  sufficiently  short,  that  the  skin  effect  may  be  appre- 
ciable. The  skin  effect  on  the  plasma  resistance  can  be  taken  into  account 
by  using  the  skin  depth  from  Table  4,  the  ratio  of  r/6  from  Table  5 and 
the  ratio  of  R/RDc  given  in  Fig.  2. 


Fig.  2 Skin  Effect  Factor 

The  resulting  resistances  for  a tungsten  plasma  cylinder  with 

three  different  radii  are  plotted  in  Fig.  3.  The  curve  line  widths 

1 8 2 n 8 

represent  variations  in  ion  density  form  10  to  10  cm  , and 
use  the  resistivity  values  for  strong  transverse  magnetic  fields. 

The  resistance  curve  for  low  temperature  is  for  the  tungsten  wire 

- - J 


r 


TABLE  4 Skin  depth. 

6 = /_n 
irfu . 

for  n.j  = 10^8  cm  8 

n.  - 1020 

T 

°K 

mm 

mm 

105 

0.560 

0.465 

106 

0.168 

0.141 

I07 

0.067 

0.058 

108 

0.015 

0.014 

TABLE  5 Ratio  of  radius  to  skin  depth,  r/S 


a 
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before  it  melts.  The  dotted  extension  of  this  curve  to  higher  temp- 

4 

eratures  is  an  estimate.  The  lowest  resistance  curve  for  the  10 

5 

to  10  temperature  range  is  for  a plasma  of  2 mm  radius  and  no  magnetic 
field.  After  the  plasma  is  compressed  and  is  in  a quasi-equilibrium 
condition,  the  kinetic  pressure  and  the  magnetic  field  pressure  are 
approximately  equal  at  the  plasma  surface.  That  is 

Ni  I? 

-7-  n + <z>]kT  = -2-7-  (7) 

nr  l 8 tt  r 

Equation  (7)  is  independent  of  r,  hence  m = 0 instability  can  easily 
develop.  If  N^/£  remains  constant,  as  it  does  in  the  one  dimensional 
code,  the  temperature  necessary  to  balance  different  currents  can  be 
calculated  from  eq.  (7)  and  are  plotted  as  the  vertical  dashed  lines 
in  Fig.  3. 

The  cross-hatched  area  on  Fig.  3 represents  the  range  of  R values 
which  would  be  obtained  if  the  plasma  reached  a semi-equi 1 ibrium  radius 
between  0.5  and  0.05  mm  with  1 to  1.5  MA  through  the  plasma.  Note  the 
1.5  n Gamble  II  impedance  is  near  the  center  of  this  range.  The  corres- 
ponding equilibrium  temperature  is  7-8  x 10^  °K. 

However,  time  integrated  X-ray  photographs  experimentally  obtained 
with  hv  > 2 keV  indicate  an  m = 0 instability  with  the  plasma  dimensions 
as  indicated  in  Fig.  4. 


Fig.  4 A portion  of  the  plasma  column  shape  as  indicated  by  X-ray  photograph 
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! The  effect  of  the  m = 0 instability  can  be  semi-quantitavely 

understood  by  considering  the  plasma  shape  shown  in  Fig.  5. 


1 


Fig.  5.  Hypothetical  plasma  column  shape 

If  there  were  no  net  transport  of  the  ions  and  electrons  in  the  z-direc- 

1 9 

tion,  the  ion  density  in  the  larger  diameter  area  would  be  ' 4 x 10 

21  -3 

and  in  the  small  diameter  area  ' 4 x 10  cm  . But  this  would  result 
in  a very  large  density  gradient  in  the  z-direction  over  a distance  of 
the  order  of  the  difference  in  radii,  so  no  transport  of  the  ions  in 
the  z-direction  is  hardly  possible  --  especially  near  the  center  of 
the  column  where  it  will  be  seen  that  the  magnetic  field  is  low.  There- 
fore, the  ion  density  in  the  small  diameter  region  is  probably  closer 
20  -3  21 

to  4 x 10  cm  than  4 x 10  . This  decrease  in  ion  density  would 

cause  an  increase  in  the  ion  density  in  the  larger  diameter  region  to 
1 9 -3 

' 7 x 10  cm  . The  change  in  ion  density  in  the  small  diameter  region 
would  call  for  an  increase  in  the  term  (1  + <Z>  )T  on  the  left  hand  side 
of  eq.  (7)  by  a factor  of  10  to  balance  the  same  current  and  its  corres- 
ponding magnetic  field.  Since  <Z>  at  the  high  temperature  will  increase 
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'25%,  the  temperature  in  the  small  diameter  region  near  the  plasma 
surface  would  be  ' 5-6  x 10^  °K.  Similarly,  the  higher  density 
in  the  larger  diameter  region  would  need  a temperature  of  ' 4 x 10^ 

°K  to  maintain  the  semi-equilibrium  condition.  These  two  conditions 
are  marked  by  points  A and  B on  Fig.  3.  If  50%  each  of  the  plasma 
column  length  is  the  small  and  the  large  diameter,  the  total  resistance 
would  be  the  average  of  the  two  values  for  A and  B,  or  approximately 
1 ohm.  For  the  observed  shape  shown  in  Fig.  4,  the  resistance  would  be 
expected  to  be  an  average  value  of  the  straight  line  through  A and  B. 

This  is  reasonably  close  to  experimental  observations. 

Also,  X-radiation  spectra  measurements  on  a plasma  similar  to  that 
discussed  here  indicate  temperatures  on  the  order  of  5 x 10^  °K  are  ob- 
tained. ^ 

For  the  conditions  along  the  line  through  A and  B the  ratio  of 
the  radius  to  the  skin  depth  6 is  ' 3.  Therefore,  current  density 
near  the  surface  of  the  plasma  is  on  the  order  of  3 times  that  at  the 
center  of  the  plasma.  This  higher  current  density,  together  with  the 
higher  resistivity  in  the  strong  magnetic  region  means  that  the  heating 
rate  near  the  surface  is  ' 10  times  that  near  the  center.  Since  the 
radiation  is  higher  for  the  hotter  regions  near  the  surface,  a temperature 
ratio  of  less  than  10  between  the  surface  and  center  regions  of  the 
plasma  would  be  expected.  Therefore,  the  ions  near  the  center  of  the 

g 

larger  diameter  regions  are  probably  less  than  10  °K. 


Since  the  plasma  at  the  surface  will  diffuse  radially  into  the 
magnetic  field,  there  will  be  a low  density  region  near  the  surface 
of  the  plasma.  Conditions  for  run-away  electrons  can  occur  in  this 
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region.  Those  that  are  produced  near  the  surface  of  the  smaller  diameter 
regions  would  move  along  the  z-direction  into  the  low  temperature  core 
regions  of  the  larger  diameter  plasma.  This  can  account  for  the  strong 
L-line  emission  which  has  been  experimentally  observed. 

A rough  check  on  the  energy  required  to  yield  the  temperatures 
estimated  above  can  be  obtained  using  the  specific  heats  given  by  Mosher. ^ 

This  is  found  to  be  ' 25  kJ,  or  a little  over  half  the  40  kJ  that  are 
deposited  in  the  plasma.  The  remainder  of  the  energy  is  radiated  both 
during  the  time  the  temperature  distribution  is  being  developed,  as  well 
as  the  time  period  after  the  semi-equilibrium  condition  is  reached. 

Since  the  above  model  appears  to  check  fairly  well  with  the  salient 
features  found  experimentally,  it  may  shed  some  light  on  the  possibility 
of  scaling  to  high  power  levels. 

Fig.  3 is  useful  to  understand  how  the  plasma  parameters  will  change 
for  increased  power  levels.  For  a current  of  10  MA  the  vertical  dotted 
line  indicates  that  the  magnetic  field  and  the  kinetic  pressure  will  be 

Q 

in  equilibrium  at  ' 2 x 10  °K,  with  no  variation  in  with  z.  With 

an  m = 0 instability  the  points  A1  and  B1  are  approximated  in  the  same 

manner  that  A and  B were  obtained  for  the  lower  current.  The  resistance 

of  the  plasma  would  then  be  an  average  of  points  along  the  A'B1  line.  It 

is  likely  that  the  equilibrium  radii  will  be  less  than  those  for  the  lower 

current,  but  it  appears  from  Fig.  3 that  the  resistance  will  still  be 

~ O.ln.  With  the  current  up  by  a factor  of  10  and  the  resistance 

2 

down  by  a factor  of  10,  the  i R heating  would  be  up  by  a factor  of  10. 

For  semi -equi 1 ibri urn,  the  radiation  must  then  be  up  by  a factor  of  10. 

j 

J 
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The  radiation  rates  are  given  byZ 

x 

Pc  = 1.62  x 10~30  /T  <Z>  ni2[<Z>  + <-T-Z--1-  q w/cm3  (8) 
and 

EZn 

Pj_  = 3.77  x 10~23  /T  <Z>  n^  < exp[-^-— ] > w/cm3  (9) 

where,  Pc  is  the  continuum  radiation,  PL  is  the  line  radiation  in  the 
absence  of  trapping  and  the  bracket  term  is  equation  8 is  given  in 
reference  7.  Since  these  terms  only  increase  as  the  /T,  and  <Z>  increase 
only  a few  per  cent,  ni  must  go  up  slightly  to  increase  the  radiation 
by  a factor  of  10.  This  is  consistant  with  the  assumption  that  the 
average  radii  will  be  slightly  smaller  than  for  the  lower  current  case. 
Therefore,  it  appears  that  if  the  current  is  increased  by  a factor  of 
10,  the  radiation  will  increase  by  approximately  10. 

However,  it  is  important  that  the  scaled-up  current  still  have  a 
time  scale  comparable  to  that  for  the  lower  current  case.  If  the  current 
rise  time  is  extended,  the  skin  depth  will  be  greater  and  the  plasma 
resistance  will  be  even  less  than  ~ 0.1  ohm.  If  the  generator  impedance 
were  reduced  by  a factor  of  10  and  the  jig  inductance  is  not  reduced, 
the  current  would  be  much  longer  in  reaching  its  value,  and  the  result- 
ing increase  in  temperature  and  radiation  would  be  less.  On  the  other 
hand,  if  the  generator  were  scaled  up  in  voltage  and  the  generator 
resistance  and  jig  inductance  held  constant,  the  current  rise  time 
would  increase  only  slightly.  Of  course,  it  may  be  difficult  to  hold 
the  jig  inductance  down  with  the  scaled  up  generator  voltage.  Also,  a 
smaller  per  centage  of  the  initial  power  supply  energy  will  be  supplied 
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to  the  plasma.  Probably  a compromise  between  a scaled-up  voltage  and  a 
scaled  down  generator  impedance  would  be  the  optimum  to  obtain  the 
higher  current.  The  determination  of  the  optimum  generator  involves 
considerable  study  which  is  beyond  the  scope  of  this  report. 


CHAPTER  II 


PLASMA  BREMSSTRAHLUNG  CALCULATIONS 

In  this  section  we  compute  the  bremsstrahl ung  spectrumfor  a high 
temperature  plasma  in  a strong  electric  field.  Our  interest  is  in 
obtaining  a theoretical  framework  for  comparison  with  the  measured 
spectrum  of  Mosher  et^  al_.6  for  an  exploded-wire  discharge.  Our  approach 
is  quite  restricted  and  is  not  an  attempt  to  describe  all  the  very 
complex  physics  of  an  exploded  wire  plasma.  Basically,  we  have  inte- 
grated the  time-dependent  electron  velocity  distribution  of  Kovrizhnykh1 0 
together  with  the  classical  bremsstrahl ung  cross  section  with  the 
results  of  Grant11  for  the  Gaunt  factor.  The  calculations  show  a 
systematic  flattening  out  of  the  spectrum  with  time  as  the  electron 
velocity  distribution  is  shifted  and  distorted  from  an  initial  max- 
wellian  shape  by  the  electric  field.  This  flatness  is  a feature  which 
is  also  seen  in  the  experimental  exploded-wire  spectrum. 

Before  describing  the  bremsstrahl ung  emitted  by  plasma  electrons 
which  have  a distribution  of  velocities,  we  consider  the  simpler  case 
of  the  bremsstrahl ung  from  mono-energetic  electrons. 

As  the  electrons  encounter  positive  ions,  they  scatter  and  emit 
bremsstrahl ung  over  a broad  frequency  spectrum.  Close  encounters 
result  in  large  scattering  angles  and  the  radiation  of  high  frequencies 
while  more  distant  encounters  give  mild  deflections  and  lower  frequency 
radiation.  To  calculate  the  radiation  one  should  use  a quantum  mechan- 
ical analysis  when  the  de  Broglie  wavelength  of  the  electron  is  compar- 
able to  the  typical  impact  parameter.  This  regime  is  specified  by 
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E > 1 3.6Z  eV, 


(10) 


where  E is  the  kinetic  energy  of  the  electron  (1/2  mv^)  and  2 the 


atomic  number  of  the  ion.  For  smaller  values  of  E a classical  analy- 


sis can  be  applied.  In  the  classical  case  the  bremsstrahlung  is  given 

11,12 


by 


X(w)  = 


«A.V 


tt/~3 


% % 01  <“„>  0)<u»t 


u=u. 


w 


G(w,v0)  (11) 


where  X(w)  is  the  cross  section  for  radiation  in  the  range  w to  w + dw. 


wb 


2 

It  has  units  of  joules  - rn  /rad. /sec.  Also  uQ  = 


. b_  = 


Ze 


H..  ^ i< 


4^omvo 


and  Hu  v,y  is  the  Hankel  function  of  the  first  kind  and  order  u_. 
o 0 


G(w,v0)  is  the  Gaunt  factor  and  contains  the  dependence  of  the  radiation 
on  frequency.  It  is  of  order  unity.  There  are  two  restrictions  on 


this  formula:  1)  w>>Wp  = plasma  frequency,  in  order  to  be  able  to  neglect 


collective  motions  of  the  electrons;  and  2)  1iw<<ffiwm=>  = l/2mv„  , since 

HiaX  0 


the  formula  is  derived  assuming  negligible  energy  is  carried  away  by 
the  photon. 

Applied  to  tungsten  wires  (Z  = 74)  for  comparison  with  the  exploded- 

wire  discharge,  equation  1 0 gives  E>76.5  keV  for  the  quantum  mechanical 

. 12  13 

regime.  Here  we  use  the  Born-El wert  formula,  ’ 


l-e‘2lTnl 


l-e“2Trn2 


In 


m(vi-vf) 

Zfiw 


2 - 


(12) 


where 


nl 


Ze 


Ze 


4^^.  * n2  = 47eo*vf  • and  vi  and  vf  are  the  electron 
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» 


speeds  before  and  after  collision,  respectively.  This  formula  allows 

2 2 

for  the  energy  carried  away  by  the  photon  (nw  = l/2mv.j  - l/2mVf  ). 

The  restrictions  are  w « wmax  and  non-relativistic  electron  energy. 

In  Fig.  6 we  have  plotted  X(w)  for  5 different  electron  energies, 
using  the  appropriate  formula,  either  equation  1 1 or  1 2 , as  indicated 
in  Table  6.  It  is  clear  from  the  figure  that  the  bremsstrahlung  spectra 
are  quite  flat  and  have  about  the  same  slope  for  all  of  the  electron 
energies  considered.  So  by  merely  observing  that  the  spectrum  is  flat 
in  a particular  experiment  one  could  not  learn  the  electron  energy. 

One  would  need  to  determine  the  high-frequency  cutoff  of  the  spectrum 
or  perhaps  examine  the  X-ray  line  emission. 

Curves  like  those  in  Fig.  6 for  target  ions  other  than  tungsten 

2 

can  be  generated  by  using  equation  11  for  electron  energies  E < 13. 62 

2 

and  equation  12for  E > 13. 62  . Since  the  equations  are  non-relativistic, 

2 

E should  be  kept  below  the  electron  rest  mass,  me  = 511  keV.  Koch  and 
13 

Motz  discuss  what  happens  if  E is  greater  than  this. 

A correction  that  has  not  been  included  in  the  spectra  of  Fig.  6 
is  that  due  to  the  screening  of  the  ion  nucleus  by  any  remaining  bound 
electrons.  We  cannot  restrict  the  analysis  to  only  fully  ionized  species, 
since  for  the  tungten  exploded-wire  plasma  this  is  not  the  case.  The 
screening  alters  the  spectra  in  the  range  of  frequencies  below  the 
frequency  w$,  where^ 


max 


zl/3 

96 


c 

v 


If  included,  it  would  merely  flatten  the  low-frequency  end  of  the  spectrum, 
whereas  in  Fig.  6 the  spectrum  gradually  rises  with  decreasing  frequency. 


-21- 


Note  that  in  the  exploded-wire  plasma,  the  electron  number  density,8 
18-1  11-1 

ng,  is  perhaps  10  cm  , so  Wp  = 5.6*10  s , and  we  easily  satisfy 
w >>  Wp. 


TABLE  6 

Data  for  Figure  6 


Electron  energy,  E 
keV 

w 

max 

sec”^ 

2 

Joule-m  -sec 

Equation 
for  Fig. 6 

5 

7.5-.1018 

1.65-10"59 

11 

10 

1.5-1019 

8.48-10"60 

11 

20 

19 

3.0*10 

4 . 37 • 1 0-60 

11 

50 

1 9 

7.5*10  y 

1 .89-1 0-60 

11 

100 

20 

1.5*10  U 

7.12-1 0-61 

12 

f 
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We  now  consider  the  bremsstrahl ung  emitted  by  a plasma  in  which 
the  electrons  have  a maxwellian  velocity  distribution.  The  bremsstrahl ung 
spectrum  is  calculated  by  integrating  equation  11  over  the  electron 
velocities,  and  the  result  is 

o 6 , -fiw/kT 

X(w,T)  = y-  y w G(w,T).  (13) 

3/3  me3  (4we0)J  kT 

Here  T is  the  temperature  and  G(w,T)  is  the  maxwell -averaged  Gaunt 

1 5 

factor  and  is  given  by  Karzas  and  Latter.  As  was  the  case  for  G(w,vQ), 
G(w,T)  is  of  order  unity.  Thus  X(w,T)  depends  on  w mainly  through 

* 

the  exponential  factor  and  so  has  a much  steeper  variation  than  in 
the  mono-energetic  case  (Fig. 6 ).  Equationl3  should  be  valid  up  to 
w : 13.6Z^/Ti,  or  Tiw  ; 74  keV  for  1 = 74. 

Applying  equationl3  to  the  tungsten  exploded-wire  plasma,  we  set 
1 = 74,  kT  = 4.4  keV  and  2.9  keV  and  obtain  the  spectra  plotted  in 
Fig.  7 . For  the  frequency  range  plotted  here  G remains  constant  at 
1.2,  so  the  spectra  are  purely  exponential.  They  are  the  same  as  the 
straight  lines  sketched  on  Mosher's  measured  spectrum. ^ 

The  objective  now  is  to  discuss  the  flat  portion  of  Mosher's  spec- 
trum which  lies  just  to  the  right  of  the  steep  thermal  portion,  above 
about  8 keV.  It  is  evidently  due  to  a stream  of  fairly  mono-energetic 
electrons  (which  we  know  yield  flat  spectra  as  in  Fig.  6),  and  these 
electrons  could  be  quite  energetic  since  the  discharge  voltage  is  nearly 
1 F}V.  We  make  use  of  the  results  of  Kovrizhnykh^0  who  gives  an  electron 


distribution  function  which  is  time  dependent  and  applies  in  the  case  of 
a strong  electric  field.  Here  we  assume  a maxwellian  distribution  at 


I 


time  t = 0, 

f<7-  t!  t , „ = WV  VZ>- 

when  the  electric  field  is  switched  on.  If  collisions  are  neglected, 
solution  of  the  Boltzmann  equation  for  t 0 gives^ 

'I-Ww?'!'  04) 

for  an  electric  field,  £,  in  the  -Z  direction.  So  the  maxwellian 
shape  moves  in  tact  out  along  the  vz  axis  at  a rate  e£/m  = Y.  The 
electrons  thus  accelerate  toward  the  anode  forming  a hot  beam,  and 
the  bremsstrahlung  emitted  as  they  encounter  ions  will  be  similar 
(although  not  identical)  to  a strictly  mono-energetic  or  cold  beam 
(Fig.  6). 

The  next  step  is  to  include  collisions.  If  the  collisions  are 
assumed  to  have  only  a small  effect,  then 

f(^’  t]  = fmax(V  VZ  ‘ Yt){  1 + 


~ fv  (VZ  - Yt)[vz  - Yt)2  + Vr2] 


1 + 3/2  [(vz  - Yt)2  - vr2] 


- 3/2  (vz  - yt)[vr2  + (vz  - yt)2]1/2 


3v  2 

-(1-  -£-)  M- 


V + Vr 


[v/  + (vz  - Tt)^]'/2  + v7  - Tt 
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► 


where  fmax(V  VZ  " = "e(|f)3/2  e"1/2s[vr  + (vZ 


yt)2] 


P = m/kT,  and  ^ = 2 10  ^2  n^2/TQ  V/ cm  with  TQ  in  electron  volts. 
This  is  Kovrizhnykh1 s result^  modified  by  us  to  apply  for  any  atomic 
number  (Z)  ion  and  any  degree  of  ionization.  Also,  we  have  let  M/m 
The  function  f(v,  t)  given  by  equation  15  is  actually  not  too  different 
from  that  in  equation  14.  Kovrizhnykh  describes  the  differences. 

The  region  of  validity  of  equation  15  depends  on  several  plasma  para- 
meters according  to  the  condition 

c »>  sk  *,(-*2X01 

m3'2^)1'2 

Unfortunately,  this  shows  that  ni  must  be  s 1014  cm-3  when  we  take, 
for  the  exploded  wire,  £,  = 105V/cm,  Z = 74,  TQ  = 4keV,  and  tmax  = 10-10 
sec-  (tmav  is  discussed  below.)  Thus  n-  is  restricted  to  values  less 

iiidA  1 

18  “3 

than  the  10  cm  reported  by  Mosher  and  so  probably  corresponds  best 
to  the  region  near  the  boundary  of  his  plasma. 

We  have  calculated  the  bremsstrahlung  spectrum  corresponding  to 
the  distribution  function  of  equation  15  . The  spectrum  is  given  by 

S(w,t)  = / f(v,  t)ni  v S (v)G(w,  v)d3v.  (16) 

Q 

S(w,  t)  has  units  of  watts/ (m  - rad. /sec). 

A sizable  amount  of  algebra  is  needed  to  put  equation  16  into  a form 
suitable  for  integration  by  computer.  See  the  Appendix  for  details. 

The  integrations  have  been  done  on  an  IBM  370  using  Gauss-laguerre 
quadrature  as  contained  in  the  library  routine  DQG  32.  The  quantity 


The  distribution  function  and  the  spectrum  depend  on  the  time, 
t,  so  that  one  must  determine  a reasonable  range  for  t.  Taking  a char- 
acteristic length  for  the  exploded-wire  plasma  as  l = 1 cm,  we  obtain 

a maximum  time,  t_„„,  from 
HldX 

Z = 1/2  (££-  ) t 2 
' 'm  ‘ max 

Thus  the  range  of  t is  0 < t < t v ~ 10~10  sec.  for  c = 105  V/cm. 

iiiaX 

When  t = t , the  electron  energy,  es£,  is  about  100  keV. 

Fig.  8 shows  f(v,t)  as  a function  of  vz  with  vr  = 0 for  t - 10~13, 
10  ^2,  and  5-10  ^sec.  We  have  taken  = 103V/cm,  I = 74, 

Tq=  4keV,  n^  = 10^  cm  3,  and  n£  = 50n^  Notice  how  f(v,t)  changes 
under  the  action  of  the  electric  field.  It  is  shifted  to  the  right 
and  slightly  distorted  near  its  peak. 

Fig.  9 shows  S(w,t)  for  t = 10"12,  10”11,  5*10-11,  and  10~10sec. 
for  n^  = 10^  ‘ and  also  10^2.  We  see  S(w,t)  flatten  out  as  t increases. 
Some  wiggles  appear  in  the  spectrum  for  the  larger  n-  - evidently  due 
to  increased  collisions.  Mosher's  measured  (time-averaged)  spectrum, 
with  both  the  steep  and  the  flat  portions,  might  be  imagined  to  arise 
from  a distribution  near  the  plasma  boundary  like  the  one  we  have 
calculated  and  a maxwell ian  distribution  deeper  inside  where  n^  is 
higher. 




w 


I 


Of  course  one  could  plot  f(v,t)  and  S(w,t)  for  other  values  of 

the  plasma  parameters  (lower  T,  for  example).  We  have  shown  only  the 

plots  of  interest  for  the  exploded-wire  plasma. 

Finally,  we  note  that  for  the  distribution  function  we  have  used, 

(equation  15)  most  of  the  electrons  may  be  classed  as  "runaway." 

The  usual  condition  for  electron  runaway  is  that  the  energy  gain  of 

the  electron  in  the  electric  field  in  one  mean  free  path  is  larger  than 

5 

its  thermal  energy.  This  gives 


’runaway 


32 

H n. 


2-4-10 
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with  Tq  in  eV.  And  for  n^  = 10  cm  , 


with  C = 10  V/ cm  we  have  Z » Z 

C - 10  < ?runaway-) 


runaway' 


Runaway  = 57  V/cm’  50  that 
(If  n.j  = 10^  though. 


APPENDIX 


Integrations  carried  out  by  computer 


On  letting  v^  = vcose,  the  electron  distribution  function  given 
by  equation  (15)  can  be  written  as  follows: 


f(v,t)  = 


ne(6/2,)3/2  e-B<a  ' b c°s9>  [1  ♦ 0 - ev2) 


+ cos  0 + (yt  - vcose)  |§  /a  - b cose 


0 2 2 2 
Bv  j.  Bv  cos  0^  e 

2 


_(1  . 


v + vcose 


J a - bcose  + vcose  - xt 


(A-l) 


where 


2 2 2 3 2 If) 

a = v + y t , b = 2vyt  and  e = 4iTe  xn  -j2^B/m. 1 


The  spectral  intensity  is  given  by 


S (w)  = / S(w,v)  f(v,t)  d v , 


where 


S(w,v)  = S G(w,v) 


cj3v  = v dv  Sineded(j> 


and 


1 1 6tt  ni2  e6 


(4TreQ)3  3/~3  m2vc3 


G(w,v)  and  S do  not  depend  on  e and  <}>.  On  carrying  out  the  angular 
integrations  we  get 


I ( v ) = /2tt  /*  f(v,t)  sineded} 


o o 
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= 2TrneC3/2Tt)3/2  [A(v)e‘(gq2)/2  + B(v)e"(f3p2)/2 

+ C(,)  r e-(ey2)/2  <iy  * ■>(*>  /“  . q)  dy 

P P 


E(v)  /q  p'ljy'VqT  **  * F(v>  /q  1^'  pHy’fq)  ^ 

P P 

(A-2) 


where 


p = V - yt,  q = V + Yt,  a = gVyt, 


A(v)  = - 


a ' 2BC  (',2  + ,Jt2  T 2vYt 


r 1 . _L 

V-  2Y2t‘ 


+ - 


seYt 


+ 7 + rir  > (' + £ ) 


bvy  t 


l 


l 


B(v)  '^^(4  + 


a 26^  v2  2Y2t2  2v^ 

+ \ ? ) £n(l  - ?r) 


«*rt  'v  Yt  ■ gvy2t2 


C(v)  - (%)  (f  - 


2.2  T T 2.2 

Y t 3vy  t 


D(v)  - (f)  tlr  - v + — V?  (a2  - to  * %)  ] 


L3v 


4vy  t 


6 ~ 62' 


E(v)  • (f)  ] 

2vy  t 


F(v)  = e/4^VY2t2 


(A-3) 


For  numerical  integration  it  is  necessary  to  modify  the  variable 
of  integration  y in  the  last  four  terms  of  I as  follows: 

Let  /e/2  y = X,  dy  = /2 7e  dx 


p 
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Also  define 

c\/m  - XU  and  p/e72  = XL. 

Then 


I =2nne(6/2rr)3/2  [A(v)e'(6q  )/2  + B(v)e"(ep2)/2 

+ m f e'*2  dx  + f Trftfmuj 


+ E(v)  fXU  x2e~X  dX  A F(v)  fXU  xV*  dX 
/I/2  TX  + XL)  (X  + XU)  (b/2)3/2  * TX  + XL)(X  + 


4 -X‘ 


XU) 


The  spectral  intensity  is  given  by 
S(w)  = / S G(w,v)  I(v)  dv 


Let  us  change  the  variable  of  integration  v to  Y as  follows: 
_ mv2  fiw 


Y = > 


2kT  " "Tf 


Then 


and 


Then 


dY  = ^ vdv 


v = /Tkf ' v fiw  , 

m [r  kT ' 


s(w)  = A0  f G[w,v( Y)]  I[v(y)]e"Y  dY 
o 


where 


(A-4) 


Ao  - 2«ne(B/2„)3/2  ^ l^)3  kT  e‘WkT 


k 


I 


Before  the  right  hand  side  of  S(w)  can  be  integrated  numerically  it 
is  necessary  to  express  the  gaunt  factor  g[w,v(Y)]  in  an  analytical 
form.  We  have  found  that  the  following  expression  for  g[w,v(Y}] 
agrees  very  closely  with  Grant's  results: 


g(w,v)  = 


rr  ln  (r78y~  ) for  y *-°08 

ln[2. 7182818  + .771588  y'-580864] 


where 


* = &T  > 

o mv 


for  y > .008 


(A-5) 


Because  of  the  presence  of  the  factor  e~T  in  egn.  (4)  the  integral 
can  be  evaluated  numerically  by  using  the  fast  and  accurate  Gauss- 
Laguerre  guadrature  scheme. 
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